Entanglement and statistics in Hong-Ou-Mandel interferometry 
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Hong-Ou-Mandel interferometry allows one to detect the presence of entanglement in two-photon 
input states. The same result holds for two-particles input states which obey to Fermionic statis- 
tics. In the latter case however anti-bouncing introduces qualitative differences in the interferometer 
response. This effect is analyzed in a Gedankenexperiment where the particles entering the interfer- 
ometer are assumed to belong to a one-parameter family of quons which continuously interpolate 
between the Bosonic and Fermionic statistics. 

Topic: Quantum measurements, Quantum optics and related technologies, Quantum entanglement and non- 
locality, New perspectives for Foundations of Quantum Mechanics from Quantum Information 



I. INTRODUCTION 



Interferometry is a widely used quantum optical technique for characterizing the non classical behaviors of light, like 
squeezing or entanglement (see for instance |l| and references therein). In particular the presence of entanglement in 
states of two multi-mode photons can be detected by means of Hong-Ou-Mandel (HOM) interferometry 0, L| • Indeed, 
the visibility of the coincidence counts at the output ports of such interferometer allow to detect entanglement in 
the incoming two-photon states: the visibility of separable two-photon inputs is in fact upper bounded by half of 
the maximum achievable visibility (e.g. the visibility associated with the entangled twin-beam state emerging from a 
monochromatic pumped non-linear crystal). In the context of mesoscopic physics Q a Fermionic analog of the 
HOM interferometer has been recently proposed in Ref. |5| : here metallic leads play the role of photonic multi-mode 
optical paths and current-current correlations play the role of the coincidence counts. Notwithstanding the change 
in the input particle statistics (from Bosonic to Fermionic) and the corresponding transition from an underlying 
bouncing behavior to an underlying anti-bouncing behavior, it turns out that the HOM interferometer retains its full 
entanglement witnessing capability. As a matter of fact the output signals of the Fermionic HOM implementation [f| 
can be easily mapped into the output signals of the standard Bosonic implementation: the only difference being 
a different definition of the "regions" which are accessible or not accessible to separable inputs. Motivated by the 
recent interest in the relations between entanglement, statistics and interferometry 0, 0, IS 0> IS & llS 111 we 
present a simple theoretical model which permits to analyze them in a broader theoretical context. In particular, 
we discuss an hypothetical HOM set-up where the two particles entering the interferometer obey to generalized 
quantum statistics [l2l Il3| . This allow us to interpolate continuously between the Bosonic HOM interferometer and 
its Fermionic counterpart. 

We start by briefly reviewing the standard Bosonic HOM interferometer. Then we introduce the quon algebra 0[l3| 
and discuss a HOM- Gedankenexperiment which operates on particles obeying to the corresponding deformed statistics. 



II. USING HONG-OU-MANDEL INTERFEROMETER TO DETECT ENTANGLEMENT IN 

MULTI-MODE TWO-PHOTON STATES 



The prototypical HOM set-up is sketched in Fig. ^ Here two multi-mode photons enter the interferometer 
following, respectively, the optical paths associated with the input ports A\ and A^. Their (possibly entangled) state 
is described by the vector 

1*) = J2 *(*i,*2)<4(fci)4(*a)|0} - (!) 

fei ,k 2 

where |0) is the vacuum and where Oj(fc) is the photonic annihilation operator associated with fc-th mode of the port 
Aj which obeys standard Bosonic commutation relations, i.e. 

[a 3 {k),a r {k')] = (2) 
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FIG. 1: Schematic of the Hong-Ou-Mandel interferometer: the two-particle state of Eq. Jy enters the set-up from the 
multi-mode input ports Ai and A2 and combines at the 50/50 beam-splitter BS. Coincidence counts JSJ are measured at the 
multi-mode output ports Bi and B2 for different values of the controllable phase-shift delay introduced in white-box area. 



In these expressions k labels the different frequencies of the modes propagating along the optical paths Ai^ (for the 
sake of simplicity the entering signals are supposed to have definite polarization) and <I>(A;i , A;2 ) is the two-photon 
spectral amplitude which satisfies the normalization condition 



£ |$(fci,fe)| 2 = l 
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Within the interferometer the modes of port A2 undergo to phase-shifts transformations introduced through control- 
lable delays (represented by the white box of Fig. [1} which transform a 2 (k) into a 2 (k)e ltpk . The two optical paths 
then interfere at the 50/50 beam-splitter BS where photo-coincidence C12 are measured at the output ports B\ and 
B 2 . We will see that the presence of entanglement in the photon input state Y$>) can be revealed by studying C\ 2 as 
a function of the controllable delays. 

Assuming perfect detection efficiency and considering the limit of long detection times T the average coincidence 
counts of the photo-detectors operating on J5i and B 2 can be expressed as 



Cl2 = lim l T dh [ T dt2 (^^ 2 (hm 
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where Ibx 2 (t) are the normalized output intensity operators defined by 

I Bj {t) E^^eH^J'^fc,)^.^) , 



(5) 
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with bj{k) being the Bosonic annihilation operator associated with the fc-th outgoing mode of port Bj. Evaluating 
the integrals and the limit of Eq. © gives 



Cia = (tf I AW I*) , 

where for j — 1, 2, Nj is the total photon number operator of the Bj output port, i.e. 



(7) 



(8) 



with rij(k) = b j j(k)bj(k) being the number operator of the fc-th mode. Analogously the average output photo-counts 
at B\ and B2 yield 



T ^°°Jo J- 



(9) 



The operators bj(k) and a,j(k) are connected through the input-output relations of the HOM interferometer, i.e. 

bi(k) = [a 1 (k) + e llpk a 2 (k)] fyft 

b 2 (k) = [ ai (k)-e^ k a 2 (k)] /y/2 , (10) 
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FIG. 2: Pictorial representation of the coincidence counts J7( dependence upon the interferometric delay tp^ in the HOM 
interferometer of Fig. Left: Bosonic case. Among the two-particle states of Eq. Q only the entangled one (represented 
by the continuous curve) can have C12 > 1/2; separable states (dotted curve) stay always below the 1/2 threshold: the gray 
area above 1/2 is not accessible to them. Right: Fermionic case. Here the area which is not accessible to separable states is 
the region below 1/2 (see Ref. Q for details). The continuous lines represent the average value l|18|l of C12. 



with (fk being the controllable phase-shift and with the coefficient l/y/2 coming from the 50/50 beam-splitter trans- 
formation. Notice that Eq. (|10|) couples only those modes of Ai and A 2 which share the same value of k and that ipk 
depends explicitly upon such index. [In Refs. 0, f° r instance it is ipk — io^t with ui^ being the frequency of the 
A;-th mode and with r being a tunable delay.] With the help of the above definitions one can show that the average 
photon number at each of the two output ports equals 1. On the other hand Eq. Q yields 

Cu = ^ , (11) 

with 

w = $*(fci,fc 2 )$(fc 2 ,fci)e l(¥ ' fc i- ¥ ' fc 2) . (12) 

It is possible to verify that for a suitable choice of the delays ifk and the spectral function $(fci, fe), the real quantity 
w can take any values in the interval [—1,1]. Correspondingly Cn can take values over the interval [0, 1]. In the 
twin-beam state case of Ref. 0, 0- for instance, the plot of w with respect to the controllable delay r shows the so 
called Mandel dip, i.e. it nullifies for r = and approaches 1 for sufficiently large r. This implies that, by defining 
the visibility V of C\i as the difference between the maximum and minimum values it assumes when varying ifk, the 
maximum achievable value of V for the two-particles state |\f r ) is 1. 

On the other hand, consider what happens when |\I>) factorizes with respect to A\ and A 2 . In this case the spectral 
amplitude has the form fc 2 ) = <j)\{ki) (^2(^2) and Eq. (|12f) becomes 

w = w sev = |^^(fc)0 2 (fc)e^| 2 ^ . (13) 
k 

Consequently, the coincidence counts C12 of factorizable input states |^) are limited into the interval [0, 1/2] and their 
maximum allowed visibility (i.e. 1/2) is only one half of the maximum achievable value (see Fig. |2J). As discussed 
in Ref. the above result generalize to any (non necessarily pure) two-photon input density matrices which are 
separable with respect to A\ and A 2 . Therefore we can use C12 as an entanglement witness for two-photon input 
states. Indeed it is sufficient to repeat the coincidence measurements for different values of the controllable delay: by 
realizing a value of C12 or a visibility V which are strictly greater than 1/2 one can conclude that the input state of 
the system was entangled. 



A. Bosons vs. Fermions 



In Ref. |5j a Fermionic implementation of the HOM interferometer has been proposed. There, effective two-electron 
states analog to the two-photon input states of Eq. (JIJ originate from a solid-state entangler y| and propagate along 



4 



metallic leads which take the place of the optical ports A\, A%, B\ and B 2 of Fig. ^ Moreover, the coincidence 
counts C12 of Eq. JSJ) is replaced by the current-current correlator of the outgoing leads Bi and B 2 . We refer the 
reader to Ref. 5] for a detailed discussion of the physical assumptions underlying such a scheme. For the purposes 
of the present manuscript it is sufficient to observe that an effective but rigorous description of such interferometer 
is obtained from the Eqs. (|TJ, Q and © by simply replacing the Bosonic commutation rules and with their 
Fermionic counterpart, i.e. 



{Oj(fc),aj/(fc')} = 
{a,,(fc),at,(fc')} = 5jj> 8 k k' 



(14) 
(15) 



with {r, s} = rs + sr being the anti-commutator of the operators r and s. With this substitution the average output 
electron-numbers of Eq. |[?JJ remain equal to 1. However the expression (fTTfl for the coincidence-counts C\i is 
replaced by 



C\2 — 



1 



(16) 



with w as in Eq. (|12|l . Since the dependence of w upon the factorizable properties of the input states is the same 
as in the Bosonic case, one can still use the visibility V as a signature of entanglement. Indeed also in the Fermionic 
case, entangled inputs can have visibility V greater than 1/2, while separable have always visibility V smaller than 
or equal to 1/2. In the case where the label k refers to the spin degree of freedom of the incoming electron this effect 
was first noted in Ref. [ll| . 

The sign difference between Eqs. (|1H and (|16fl implies that separable states of Fermions are forced to have C12 
greater than or equal to 1/2 while separable states of Bosons have C12 smaller than or equal to 1/2 (see Fig.|2J. This 
phenomenon can be interpreted in terms of the different bouncing and anti-bouncing attitudes of Bose and Fermi 
statistics. Indeed due to the exclusion principle one expects the coincidence counts C\i of Fermions to be typically 
higher than the corresponding Bosonic coincidence counts. To make this a quantitative statement let us consider the 
average value of Cn over all possible two-particle states of Eq. Q), i- e - over & U the normalized two-particles spectral 
amplitudes k%). A simple symmetric argument can be used to show 



w = I dfi($) 



J2 $(fei,/c 2 )$*(fe,fci)e i( ^-^) 



l/M , 



(17) 



with M being the number of the fc-th modes. Therefore the average coincidence counts is 

l/2-l/(2M) Bosons 
l/2 + l/(2M) Fermions. 



C12 



(18) 



It is worth noticing that entangled input states are, to some extent, insensitive to the complementarity behavior of 
Eq. I|18[) : indeed, as shown in Fig. [3 there are no regions which are strictly forbidden to them by the particle statistics. 



III. HOM INTERFEROMETRY WITH QUONS 



In this section a Gedankenexperiment is analyzed where the particles entering the HOM interferometer are assumed 
to be quons |l2t IT^ | . These have been introduced by Greenberg as an example of continuous interpolation between 
Bose and Fermi algebras (see Refs. 0,^3 f° r alternative definitions). We will use them to describe the sharp transition 
from the left part to the right part of Fig. [21 in terms of a continuous deformation of the particle statistics. 

For q E [—1,1] the quon algebra is obtained by replacing the relations and (|15fl with the identity 

Oj (fc)at, {k')-q at, {k')a 3 (k) = 5 rj , 5 kk , , (19) 

and by neglecting the relations @ and (|14[1 . For \q\ < 1, Eq. (|19fl can be used to define a valid non-relativistic 
field theory but poses serious problems in deriving a reasonable relativistic quantum field theory |l2l . A proper 
Fock-like representation can be derived upon enforcing the vacuum condition 



a i (fc)|0)=O, 



(20) 
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FIG. 3: Plot of the maximum and minimum values of C12 of Eq. i'2'M allowed for a two-particle quon state as a function of the 
parameter q of the quon algebra. The shaded region represents the allowed values 1241 of C12. The gray region is accessible 
only to non factorizable input configurations - see Eq. J2SJ- The continuous line represents the average value l|2(il of C12 over 
all possible input two-particles states |^). F or q — 1 and q = — 1 the boundaries coincides with those of Bosonic and Fermionic 
case respectively. The Boltzmann statistics |l2| q = admits only the value 1/2. 



for all j and k |l2l Ha |. In particular, it can be shown that for \q\ < 1 the squared norms of all vectors made by 
polynomials of at(fe) acting on the vacuum state |0) are strictly positive. For q = 1 and q = — 1 instead the squared 
norms are never negative and nullify for those configurations which are, respectively, totally symmetric and totally 
antisymmetric under permutations. This allows us to recover the Bosonic and Fermionic statistics as limiting cases 
of Eq. (|19|l without explicitly imposing the conditions J5J and flj 1 4|> . respectively. Of particular interest is also the 
q = case whose statistics can be interpret jl2j as quantum version of the Boltzmann statistics, based on a system 
of identical particles having infinite number of internal degree of freedom. 

The possibility of defining a proper Fock-like structure for the quon algebra l|19[) implies the existence of number 
operators n,j(k) which satisfy the standard commutation relations 

[rij(k),aj>(k')] = -6jj> 6kk> aj(k) , (21) 

and which reduce to a^(k)aj(k) in the Bosonic and Fermionic limits 0, ^1- In the following we will adopt a 
pragmatic point of view assuming that the equations of the previous section yield a legitimate description of our 
quon HOM interferometer. This is in part justified by the fact that the input-output relations (|10|) map the quon 
annihilation operators a,j(k) into annihilation operators bj(k) which still satisfy the quon algebra l|19[) and the vacuum 
condition (|20|l . The only technical problem of our approach comes from the fact that for \q\ < 1 Eq. (JJJ is not 
necessarily a real quantity (it is a consequence of the fact that no commutation relation between a\{k\) and 02(^2) is 
defined [lql. Consequently we replace the product N1N2 with its Hermitian part, i.e. 

C12 = (*| 1 - 2 2 2 1 |tt) = Sfte (*\N X N 2 \*) . (22) 

[In this expression jVi^ is defined through Eq. JSJ in terms of the number operators n-j(k) of Eq. J2J]. The right 
hand side term of Eq. (1221) can now be computed by inverting the transformation IjlUI) and expressing the input state 
|\&) of Eq. in terms of the operators bj(k). With the help of the relations (|21|) and (|19H and using the vacuum 
condition (|20|) one can then verify that the quon "coincidence counts" of the g-algebra are [lq . 

CM = ^P (23) 

with w as in Eq. (|12[). Analogously one can verify that for normalized input states Eq. J3J still holds and that the 
average number ix,2 is equal to 1 for all q. As expected, for q = ±1 Eq. (|23|l reduces to the Bosonic and Fermionic case. 
More interestingly for q = 0, Cuiq) does not depend upon the two-particle input state |\1>) and has constant value 
1/2. This is exactly what one would expect from a classical model of the interferometer confirming the Boltzmann 
interpretation of the q = algebra. 
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Taking into account that for generic input state one has w e [—1,1] the following bounds can be derived: 

< Cm(q) < ^ forge [0,1] 

(24) 

^<C 12 (g)< V forge [-1,0]. 

On the other hand for factorizable amplitudes <&(fci,fc2), the function iu obeys Eq. Therefore for these states 

one has 

^<C 12 (g)<| forge [0,1] 

(25) 

± ^C 12 (g) ^ ±^ forge [-1,0]. 

Finally, we can use Eq. (|17|l to compute the average value of value of C12 with respect all possible two-particle input 
states \$f), i.e. 

C i2 (q) = l/2-g/(2M) . (26) 

The above constraints and Eq. I|26|) have been plotted in Fig. |3J they indicate that the transition from the Bosonic 
regime to the Fermionic regime is characterized by a critical point at g = 0. Here the values C12 > 1/2 which for g > 
were pertinent to non factorizable input states |\&), become accessible to factorizable configurations. At the same 
time, however the values C\i < 1/2 become inaccessible to them. Following the discussion of the previous section 
this effect can be interpreted as a continuous transition from a bouncing behavior to an anti-bouncing behavior with 
q = corresponding to the classical "neutral" point. Moreover, Eq. I|2t)|) provides an indirect confirmation of the 
interpretation of the q = case in terms of the statistics a system of identical particles having infinite number 
of internal degree of freedom. In effect according to such expression the limit M — ► 00 and q — > of the average 
coincidence counts Ci2(g) are equivalent. 



IV. CONCLUSIONS 



We analyzed how different quantum statistics affect HOM interferometry. In particular we focused on the in- 
terferometer response to initially separable/entangled inputs. By introducing a family of g-dcformcd algebras we 
interpolated between the Bosonic and Fermionic regimes which were previously discussed in Ref. [BJ . Our results indi- 
cate a progressive attenuation of the bouncing behavior when moving from the Bose statistics to the Fermi statistics 
a progressive attenuation of the bouncing behavior. 
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